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Abstract 

We show a new class of embedding solutions of D5 brane, which wraps on S 5 in 
the AdSs x S 5 space-time and contains fundamental strings as U(l) flux to form a 
baryon vertex. The new solution given here is different from the baryon vertex since 
it consists of two same side (north or south) poles of S 5 as cusps, which are put on 
different points in our three dimensional space. This implies that the same magnitude 
of electric displacement exists at each cusp, but their orientations are opposite due 
to the flux number conservation. This configuration is therefore regarded as a D5-.D5 
bound state, and we propose this as the vertex of a baryonium state, which is made of 
a baryon and an anti-baryon. By attaching quarks and anti-quarks to the two cusps 
of this vertex, it is possible to construct a realistic baryonium. 
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1 Introduction 



In the context of string/gauge theory correspondence [H El E], the baryon has been 
studied as a system of fundamental strings (F-strings) and D5-branes wrapped on S 5 
in AdS 5 x S 5 space-time gQEl El El El HH H21 IEI - They correspond to quarks and the 
baryon vertex respectively. The F-strings are partially dissolved as a U(l) flux in the 
D5 brane, and their remaining parts flow out from one (or two) cusp(s) on the surface of 
the D5 brane as separated free strings. The baryon vertex has complicated structures 
which are given as solutions of the equations of motion for the D5 brane embedded in 
an appropriate background, which is dual to the confining gauge theory (for example 
[T5l [T6| [T7]). This picture has been recently studied furthermore [18] along the Born- 
Infeld approach given in [5]-[IT]. and also extended to finite temperature theory [T9|l20] 

Here, we show new kinds of configurations, which are obtained as solutions of the 
same equations with the one which gives baryon vertex solutions. But the new solutions 
given here are different from the baryon vertex. The baryon vertex wraps whole S 5 
once. Namely, it covers all range of the polar angle (9) of the S 5 , < 9 < ir, once. 

On the other hand, the configuration of the new solution covers twice one polar side, 
for example in the range of 9q < 9 < tt where < 9 . And, at any 9 in this range, this 
configuration exists at two different point in our three dimensional space. Then this 
configuration looks like a line with a finite length. On this line, the point of 9 = 9q is 
at its center and the two end points are given by the same 9 = tt. In other words, it 
starts from one end at 9 = tt and arrives at 9 = 9$ (the center of the configuration) 
then goes to the other end point 9 = tt along another half path. This solution can 
be interpreted as the connection of two U(l) fluxes with opposite charge. This fact 
implies that it describes the same polar sides {9 = tt) of connected two S 5 s wrapped 
by D5 and anti-D5 (-D5) branes respectively. 

This can be regarded as a D5/D5 bound states. Similar D-brane embeddings have 
been found for different D-branes in different backgrounds [2H [221 E3]. The essential 
point of such solutions is that the configuration covers two different points in our space 
at the same point of a world volume coordinate (here 9) of the D-brane. In this sense, 
our solution is essentially the same type with the former examples. 

The two end points are the cusps, where opposite sign of U{1) fluxes exist. Then the 
F-strings attached at these two cusps have also opposite orientations with the same 
number. This is considered as the baryonium or the bound state of a baryon and an 
anti-baryon. 

The energy and the configuration of this baryonium vertex depend on boundary 
conditions of the equations of motion. So, varying the boundary conditions, the relation 
between the vertex energy and the distance of the two cusps is examined. And we could 
find a minimum vertex energy at a finite distance between the two cusps. This implies 
that the baryon and anti-baryon bound state is stable against vanishing to the vacuum. 
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In Section[2]we give our model and D5-brane action with non-trivial U(l) gauge field. 
And the equations of motion for D5 branes are given. In section El we give D5/D5 
solutions as baryonium vertex. And its configuration and energy, which depends on the 
configurations, are examined. In the section 4, the differences between the baryonium 
and split baryons are discussed. And in the final section, we summarize our results 
and discuss related directions. 



2 Model 

2.1 Bulk background 

We derive Dh-Dh solutions as baryonium from the equations of motion given by the 
action of D5-brane which is embedded in a supersymmetric lOd background of type IIB 
theroy. The background solution should be dual to the confining gauge theory since 
the baryonium examined here is a bound state of quarks. While there may be some 
such solutions, we consider the following background [T5l [16l [T7] . 

ds\ = e*/ 2 [^ r)ilv dx^dx v + ^dr 2 + R 2 dn 2 A , (1) 

which is written in string frame. At the same time, the dilaton $ and the axion x are 
given as 

e* = l + ^, X = -e-* + X o, (2) 
and with self-dual Ramond-Ramond field strength 

G(5) = dC {4) = AR A |vol(S 5 )d#i A ... A d6 5 — ^dt A ... A dx 3 A dr^j , (3) 

where vo^S* 5 ) = sin 4 9ivo\(S 4 ) = sin 4 9\ sin 3 Q 2 sin 2 #3 sin # 4 . 

This solution, ([I])-©, is useful since the confinement of quarks are realized due to 
the gauge condensate q = (F 2 U ) [16],[T7j, which is given by the coefficient of 1/r 4 for the 
asymptotic expansion of the dilaton at large r. And furthermore, M=2 supersymmetry 
is preseved in spite of the non-trivial dilaton is introduced. We can assure through the 
Wilson loop that q 1 ^ 2 is proportional to the tention of the linear rising potential between 
the quark and anti-quark [T7]. In the present case, q is essential to fix the size of the 
baryonium and stabilize it energetically as shwon below. 

We notice that the axion \ corresponds to the souce of D(-l) brane and it is Wick 
rotated in the supergravity action. This is necessary to preserve the supersymmetry. 

2.2 D5 brane action 

The baryon is constructed from the vertex and N fundamental strings, and the vertex 
is given by the D5 brane wrapped on the S 5 of the above metric. The iV fundamental 
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strings terminate on this vertex and they are dissolved in it [H [5] as U(l) flux. The 
D5-brane action is thus written as by the Dirac-Born-Infeld (DBI) plus WZW term [7] 



S D5 = ~T 5 J d^e^^-det (g ab + F ab )+T 5 J d 6 £i (1) A C {5) , (4) 
gab = 9 a X di,X Gfj, u , C 0l ... O5 = (9 ai X M1 . . . d a5 X fl5 G fll ,,, fl5 . 



where F a b = 27ra'F a b and T 5 = l/(g s (27r) l s ) is the brane tension. 

The D5 brane is embedded in the world volume £ a = (t, 9, #2, ... , 9 5 ), where (# 2 , . . . , # 5 ) 
are the S 4 part with the volume of f2 4 = 87r 2 /3, where we set as 9 1 = 9. Restrict our 
attention to SO(5) symmetric configurations of the form r(9), x(9), and A t {9) (with 
all other fields set to zero). Then the above action is written as 

S = T 5 fi 4J R 4 J dt d9{- sin 4 9 ^Je* (r 2 + r' 2 + (r/R) V 2 ) - F t 2 d - F te D}, (5) 

where the WZW term is rewritten by partial integration with respect to 9, and = 
87r 2 /3 is the volume of the unit four-sphere. The factor D{9) is defined by 

d e D = -4 sin 4 9 , (6) 

and is related to F t g by the equation of motion for A t as 

D — , (7) 



We call this D as displacement, and it is given by solving ([6]) as follows , 
D(u,9) 



3 3 

— (vk — 9) ^ — sin 9 cos 9 + sin 3 # cos 9 
2 K J 2 



The meaning of the integration constant, defined in the range of < v < 1, is given 
below. 

Next, the action is rewritten by eliminating the gauge field in terms of (GO) to obtain 
an energy functional of the embedding coordinate only 0: 

U = -4-: / dd e^^Jr 2 + r' 2 + (r/R) V 2 JvJ&) . (9) 

^(0) = D(u, 9f + sin 8 (10) 

where we used T5Q4R 4 = N/(3Tr 2 a'). Using this expression ([9]), we consider the mean- 
ing of the integration constant v given in ([6]). In the below, we solve the equation of 
motion for r(9) and we find that it has two cusps or singular points at r{9) = r(ix) and 



1 U is obtained by a Legendre transformation of L, which is defined as S — J dtL, as U 

L 

dF tt 



jrk-Fte — L. Then equations of motion of (0) provides the same solutions of the one of U . 



r(0), namely at 8 = it and 9 = 0. At these points, r' = dgr diverges and x' ~ for 
q = 0. The configuration near these positions represents the bundle of the fundamental 
strings. The numbers of the fundamental strings at the cusps are estimated as follows. 
At 9 = tt and for q = ($ = 0), we obtain the following approximate formula 



Since J dr represents the bundle of a fundamental string, the total number of 
fundamental strings is given by N, which are separated to iV(l — v) and Nv to each 
cusp point. The meaning of v is then the ratio of this separation, so it must be defined 
as < u(= k/N) < 1, where k(< N) is an integer. 

By the definition of U, Eq.([H]), U is positive, and it is proportional to \D(v, 0)| or 
\D(u, 7i") | . Then the total number of the flux is counted as N when we sum up the one 
of the two cusps at 9 = and 9 = tt. However, we notice here the orientation of the 
flux of U(l) current, then D defined by (JZ|) could takes two possible value, D = ±\D\, 
depending on the orientation of the flux. For the case of opposite orientation, the total 
flux number would be counted as — N. This is regarded as the anti-baryon vertex. 

Then two possible flux numbers are assinged as ±iV(l — v) and ±i>N at each cusp. 
For the split baryon, which extends between the cusps at 9 = and tt, we find ±iV 
since the baryon must be a color singlet. However, we found new solutions, which 
extend between the cusps at the same 9(= or tt) as shown below. In this case, for the 
solution with two cusps at 9 = n, we must choose the flux-combination as ±N(1 — v) 
and ^N(l — v). And for the one with the cusps at 9 = 0, the flux should be assigned as 
±uN and ^fuN. Then the total flux is zero in both cases, and we call these solutions 
as baryonium. 

In the case of q > 0, r'(7r) or r'(0) does not diverges any more, but these points, 
9 = tt and 9 = 0, are singular and tension to deform the D5 brane is observed. This 
tension can be cancelled out by adding the fundamental strings whose number is given 
by the bundles observed for q = 0. This picture is very naturally understood since the 
system of the brane and the fundamental strings deforms continuously with the scale 
parameter q. 

Indeed above mentioned statements can be explisitly checked. For this purpose we 
calculate the tensions at the cusps [HI [191 EZ] in Appendix B. To compare tensions, we 
take a "vertical limit", namely, r[ — > oo, — > oo. Then the following equality holds; 




(11) 



And similary, we obtain the following at 9 = 0, 




(12) 



6U_ 

Sri 



(1 - u)N 



5U F 



5ri 



(13) 
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The tensions in the ^-direction vanish in the vertical limit. The above equality means 
the tension of the cusp equals to (1 — v)N times that of F-string automatically in the 
vertical limit in the case of q > 0. Similary, we find the vN times tension of F-string 
on the other cusp. 

3 Baryonium States 

Equations of motion 

In terms of (Q, we could obtain two kinds of baryon configurations [18]. In 
both cases, we should notice that the solutions r{9) and x{9) cover whole region of 9, 
< 9 < 7r. However there are other kinds of solutions, which cover only a part of the 
variable 0, i.e. (i) 9 < 9 < tt or (ii) < 9 < 9 X , where 9 > 9 X . 

Here we identify these type of solutions embedded in (i) or (ii) as the baryonium. 
Namely, for the baryonium solution, 9 does not cover all the region for both cases. 
For example, in the case of (i), the solution x(9) starts at x(7t)(= x_ < 0) and passes 
x{9q) = smoothly, then arrived at x{ti){= x + > 0). So, this configuration extends 
from X- to x + in our real three space. Here we should notice that the two end points 
are at 9 = n then the flux numbers at these points must be the same but with opposite 
direction. This corresponds exactly to the baryonium as mentioned in the introduction. 
Similarly, we can consider the baryonium of region (ii), whose end points are at 9 = 
with different flux numbers. 

The difference of the baryon and baryonium solutions is reduced to the difference 
of their boundary conditions. Generally speaking, in solving the differential equations, 
the boundary condition determines the integration constants which correspond to the 
constant of motion like energy. In the present case also, we introduce such a constant, 
which discriminates the solution of baryon and baryonium. 

In order to introduce such a constant, we rewrite (Q by changing the integration 
variable from 9 to x as follows 

U = -^L f dx e^ l2 \jr 2 9 2 + r 2 + (r/R) A sjvjff), (14) 

where dots denotes the derivative with respect to x. We can introduce an integral 
constant h as a "Hamiltonian" for the corresponding "time variable" x as follows 

h = rp r + 9p e -L, (15) 

where 

L = e* /2 ^fr 2 9 2 + f 2 + (r/Ry sjv v {9) (16) 



5 



dL n 
Pr = -wr = rQ 

or 



86 



and 



\ 



[IT) 



(18) 



Then h is written in terms of the momentum as 

and the equations of motion are obtained as 

• P r a : 
r = — , a 



+ P 2 r , 



,2 ^ r 



r 2 Q ' 



(19) 



(20) 



(21) 



These equations are convenient to find the baryonium vertex solution as seen below. 

20- 




Fig. 1: V u {0) for v = 0.2, 0.4, and 0.5. The horizontal line shows a sample line of V v { 
= 2.0. The crossing points between this line and the curve of V v represnets 6q. 



Before giving explicit solution, we show how the type of solutions is controlled by h. 
From the definition of h given in Eq. (j!5p . we obtain 



then from the reality of the solution, the next constraint is obtained 

r4 + <? T , ^ , 2 
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V v > h 2 . (23) 



At the end points, x = x±, 9 = it and r could become large, however r takes its 
minimum r ~ ro at the mid point x — 0. Expressing as 8\ x= q = 9q at this point, the 
value of 6*o is obtained by solving the equation, 

^K,, = h> . (24) 

This equation has two solutions for 9q when h is large as shown in the Fig. [TJ but 
there is one or no solution for small h. The solutions of the former (latter) case of large 
(small) h are identified with the baryonium (baryon). 



Consider the baryonium solution. From ff24|) . r is given as 



V, 




r o = T777TT -1) • ( 25 ) 



While the exact solutions are obtained by solving the above four Hamiton equations 
(1201) and ( T2~TT) as given below, we show an approximate solution in the region where we 
can assume r ~ and r ~ ro near x = 0. Under this assumption, we obtain from (1221) 



= ±^yfiUqy/Vu(B)-V v (6 ) (26) 

then we solve this as 

x{6) = ± / dO . (27) 

ha r y/ro 4 + qy/V v (0) - V v %) 

This solution is symmetric with respect to x = axis in x-9 plane. The important 
point of this approximate solution is that the solution runs from x = to two opposite 
directions, however they are going to the same pole on S 5 but with different values of 
x. In order to see the behavior of the solution far from 8 , we must solve the exact 
form of equations. Actually we can find the solutions as stated above exactly, namely 
they satisfy this symmetry at all 9 even if the assumption imposed here is not satisfied. 

Which pole, 9 = or n, is chosen depends on the value of 9q. Notice that V u {9) has 
a minimum at 9 C which is given as a slolution of 

7tv = 9 C — sin 9 C cos 9 C . (28) 

and we find the minimum value as V U (9 C ) = sin 6 (6 l c ). This implies that the pole 9 = 
(9 = 7r) is chosen for 9 < 9 C (9 > 9 C ). This is understood well from the Fig. [TJ We 
notice however that the number of 9q depends on the value of h as seen for v = 0.2 as 
seen in the Fig. HJ The situation is however changed by the value of h to find two 9q 
also for v = 0.2. 
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But for the case of v = 0, there is only one 0o f° r an Y value of h. In this case, the 
baryonium is constructed by iV-quarks and anti-iV-quarks attached at each end points. 
In other words, we obtain a bound state of baryon and anti-baryon without any loss 
of quark and anti-quark by their pair annihilation. 

On the other hand, for the case of v = 1/N,we expect another interesting baryonium 
configuration which is constructed by one quark and one anti-quark. This state is very 
similar to the usual mesons, but it is different from them in the point that the D5 
vertex is included except for a quark and an anti-quark pair in this state. 

Numerical solutions 

Here we give the explicit baryonium solutions mentioned above. The equations are 
solved numerically since it would be impossible to solve them analytically. Firstly, we 
give a way to obtain the symmetric solutions. It is convenient to set the following 
boundary conditions at x = as 

0(0) = O , r(0)=r , and p r (0) = p*(0) = (29) 

They are given as follows. First, an appropriate 0o is given for fixed h, then ro is 
determined by the above relation ([25]) . The last two conditions, p r (0) = pe(0) = 0, are 
necessary to obtain a baryonium solution which is symmetric with respect to x = 
axis. 



r(x) 




-1 -0 . 75-0.5-0 .25 0.25 0.5 0.75 1 ' 



.75-0.5-0.25 0.25 0.5 0.75 



Fig. 2: The typical D5/D5 solution for q = 0.3, u = 0.5 at h = -1 and 6 = 1.6708. The 
boundary conditions are p r (0) = po(0) = 0. 



An explicit example of the baryonium solution is obtained for v = 0.5, h — — 1 and 
0o = 1.6708 under the above boundary conditions (T29]) . The results are shown in the 
Figs. [2] and [3j From the Fig. [21 we can easily understand the fact that this solution is 
interpreted as the D5 / D5 bound state solution. This kind of solutions can be obtained 
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Fig. 3: The 3d graphic of D5/D5 solution given in the Fig. [2]for q = 0.3, v = 0.5 at h = — 1 
and 6»o = 1.6708. 



at any u, and it can be considered as the vertex part of a baryonium of (1 — v)N quarks 
and (l — i/)N ant-quarks which are attached at each end point of the vertex for general 

/A 

Here, in the Fig. HJ we show the u = 0.7 case as an example. 




-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 -0.6 -0.4 -0.2 0.2 0.4 0.6 



Fig. 4: D5/D5 solution for q = 0.3, v = 0.7 at h = — 1 and 9q = 1.9962. The boundary 
conditions are p r (0) = pe(0) = 0. 



Asymmetric solution 

While the symmetric solutions are examined above, the various asymmetric solutions 
are also obtained when we solve the same equations with a slightly different boundary 
conditions from the one of the symmetric solutions. For example, they are obtained 
by changing the boundary conditions, ( !29l . These solutions are also regarded as the 
baryonium since they connect two cusps in the same side, at 9 = n (or at 0) at different 
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X. 

An example of a little asymmetric solution is shown in the Fig. [5], where boundary 
conditions are changed to p r (0) = and p$(0) = 10~ 4 . 



r(x) 



12 




1.5 



1.0 





-0.5 



0.0 



0.5 



-0.5 



0.5 



Fig. 5: The typical D5/D5 solution for q = 0.3, u = 0.5 at h = -1 and 6 = 1.6708. The 
boundary conditions are p r (0) = 0, pe{0) = 10 -4 . 

In general, the energy of the asymmetric solution becomes higher than the symmet- 
ric one. In this sense, the stable configurations can be considered as the maximally 
symmetric one. Then, we consider hereafter symmetric solutions. 

Another kind of solution obtained by a different boundary conditions is the one 
called as the split baryon, which represents a baryon vertex. We give comments to this 
solution in the next section. 

Stability of the baryonium vertex 

In the next, we study the stability of the baryonium solution obtained above. From 
the viewpoint of energy, we concentrate on the symmetric vertex solutions, which are 
discriminated by its length L. Here L is defined as 



where we assume x + (tt) > x-(tv). Depending on the boundary value 6q, both the L 
and the vertex energy U, which is given by (}Tlj) . vary. So, by varying 9 for fixed h, 
the relation between U and L is examined for the symmetric solutions. This relation 
could give us a critical check for the stability of the baryonium configuration. 

The numerical results of U and L for h — — 1 and v = 0.5 are shown in the left hand 
side of Fig. El where the normalized energy density, U n = U / ( 3 ^ a , ) , is shown instead 
of U. From this we can see that a clear minimum is found at definite and finite value 
of L. Namely, the state for large L and also for small L needs large energy to form it. 
Especially, the latter fact means that the baryonium vertex does not vanish to L = 0. 
For other values of h, similar U-L relations are expected and this is assured as follows. 




(30) 
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Fig. 6: The left: The normalized vetex energy U n = U /( 37T -2 a > ) versus its length in the x- 
direction for D5/D5 solution for v = 0.5, q = 0.3, R = 1 and h = — 1. The line represents 
U n = 1.48L. The right: The equi-potential curves in (—h,L) plane. The point A denotes the 
bottom of the potential. 



Here, U depends on the two parameters h (or tq) and 6q except for the external 
parameters R and q. Then we can write as U = U(h,r ) or U = U(h,6o). Meanwhile 
it is possible to replace one of the parameters by a physical quantity, for example L, 
as U= U(h,L). 

In order to assure the minimum of U in the two dimensional parameter space, we 
show the equi-potential curves (contour for U) in the h-L plane. The numerical results 
are shown in the right hand side of the Fig. [6j From this, we find a minimum does 
exist near the point A, and we find that the minimum given in the left figure is near 
the same L of point A. In general, we can consider many kinds of paths in this plane 
to study similar U-L relation along it. 

From the above results, we could read the following points. 

(i) The energy density U has a minimum at a finite L. This implies the stability 
of the baryonium configuration against the tachyon, which would appear as strings 
connecting between D5 and anti-D5 branes. The mass square of this string is given 
as L 2 /R 4 — l/(2na') [23], then it is stable for enough large L. Here the size of L is 
measured by the scale R, the AdSs radius, so we can set the parameters to satisfy the 
inequality, 
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at least around the L where U is minimum. 

(ii) Secondly, at large L, U increase linearly with respect to L, and we can approxi- 
mate its behavior as 

U = r B L. (32) 

where tb denotes the tension of the baryonium vertex. The value in the case given in 
the Fig. [6] is evaluated and shown by a line. However, we must notice that the above 
result is obtained for h = — 1 and t# depends on h as seen from the right hand figure 
of U(h,L). So we must be careful to study the tension of the baryonium vertex. On 
this point, we do not discuss furthermore. 



4 Baryonium and split baryon 

As discussed in [HI HE], the baryon vertex has two types of solutions, that is, the point 
and split one. The latter is similar to baryonium vertex in the sense that it has a finite 
length in our space. 

As mentioned above, the baryon vertex configurations are also obtained by solving 
the same equations with somewhat different boundary conditions. As given in [TBI [19] . 
an easy way to obtain the baryon vertex, which extends also in the x-direction with 
the two end points x(n) and x(0), is as follows. First, set the boundary condition at 
9 = 9 C , where 9 C is the minimum point of V v {9) as given in ( 1281) . as 

Pe\e=e c = V ^ , p r \e=e c = . (33) 

This boundary condition is necessary to embed the D5 brane on the whole region of 
S 5 . In other words, the polar angle 9 must cover the whole range < 9 < tt, and 
this becomes possible under the condition rj ^ in fl33|) . This is the important point 
to identify the obtained configuration with the baryon vertex. Namely, the two end 
points are at different polar angles, 9 = and 9 = n. 

For the split baryon vertex obtained in this way, we examine its energy. Roughly 
speaking, its configuration can be separated to the two parts of extending to x and to r 
directions. And each part shears the energy of the vertex. We find that the split baryon 
vertex can be smoothly deformed to the point vertex by suppressing its energy. Then 
the configuration of minimum energy is realized by the point vertex solution r = g 1 / 4 
[TBI [T9] . This point is the important difference from the case of the baryonium, which 
could not be pushed to a point in spite of the fact that the baryonium is made of 
D5/D5. 

As for the total mass of the baryon or baryonium, we must add F-strings under the 
appropriate conditions called as no force condition [18]. In this case, we could find 
that the minimum energy of split baryon is realized when the length of the F-strings 
vanishes but the vertex length is finite. 
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Similar situation is also expected for the case of the baryonium. We examined the 
energy of baryonium for v = 0.5, q = 0.3, R — 1, r max = 10 and h = —1. Then we 
could find that the minimum of the energy is found when the F-string length vanishes 
as in the split baryon. However the details of the analyses are not given here. We will 
give them in the near future. 

5 Summary and Discussion 

We find a baryonium solution, which can be interpreted as bound state of baryon and 
anti-baryon, by solving the equations of motion for the D5 brane action. The reason 
why the bound state D5/D5 is obtained from the D5 brane action is that the action 
used to solve the equations contains the displacement flux operator D in the squared 
form D 2 . This fact enables us to obtain the baryonium configuration, which is made 
by connecting a D5 brane and its anti-brane, from a D5 brane action which we used. 

The configuration of the baryonium vertex looks like a string in our three dimensional 
space and we find that its size or length L. Its energy U depends on the length L, and 
we can show the minimum of U is realized at finite L. Then the size is kept finite in its 
stable state, and we could assure that its configuration in the bulk is in a form that the 
D5 and D5 are separated enough not to be destabilized by the tachyon. Then we can 
say that the baryonium state found here would be stable and it would be difficult to 
observe its decay into mesons. This would be related to the selection rule for hadronic 
decay and the resultant narrow width of the baryonium |25j. 

The baryonium vertex solution given here is similar to a baryon vertex configura- 
tion which also looks like a string in our three dimensional space. And they both 
are obtained from the same equations of motion. However, one of the values of the 
displacement D at the two cusps for the baryonium is different from the baryon. The 
sum of D at the two cusps is zero for the baryonium, but it is iV for the baryon. The 
energy minimum for the baryon is realized for L = 0, which looks like a point in our 
space. Meanwhile, its configuration looks string like in the bulk. 

A real baryonium should be made of the vertex and fundamental strings attached at 
cusps as performed in the case of baryons. Then the mass spectra for the baryonium 
are examined to compare the spectra given in recent experiments [26] . The tetra-quark 
meson corresponds to the baryonium of v = 1/3 and N = 3. We can estimate the mass 
spectrum of this state. This would be given in the near future. 

Appendix A 

Here we show another formulation of solving the equations of motion derived from ([9]) 
used in [Tij [18]. Firstly, rewrite (Q in terms of a general worldvolume parameter s 
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defined by functions 9 = 9(s), r = r(s), x = x(s) as: 

N 



U 



3ir 2 a 



- J ds e* /2 y/r*9 2 + f 2 + (r/R)*x 2 ^Jv u {9), (34) 



where dots denote derivatives with respect to s. Then the momenta conjugate to r, 9 
and x are given as 



Pr = rA, Pe = r 2 9A, p x = (r/R) 4 xA, A = e* /2 _ 

'r 2 9 2 + f 2 + (r/R) 4 x 2 



Since the Hamiltonian that follows from the action f)T4j) vanishes identically due to 
reparametrization invariance in s. Then we consider the following identity 



2H = p). + V \ + ^p 2 x - {V v {9)) e* = . (36) 



Regarding this constraint as a new Hamiltonian, we obtain the following canonical 
equations of motion, 

r = p r , p r = ^p 2 x R 4 + V f,+\ {V v {9)) e*d r $, (37) 

9 = p e = -6 sin 4 (vrz/ - 9 + sin9cos9) e* (38) 

* = ~TP*> = (39) 

The initial conditions should be chosen such that H = 0. By solving these equations, 
we could find the same solutions given above. 



Appendix B 

Here we calculate the tensions at the cusps [TSJ [THl EZ] • In the present model the tension 
has r-component and x-component generally. Denoting T\ = r(ir) and x\± = x±(iv), 
where Xi_ < 0, and x 1+ > 0, the tensions are given by, 



SU {l-v)N e^' 2 r' x 5U {l-v)N e^' 2 {r x / R) A x' 



i± 



1± 



Sn 2na' y/rf + r? + {r 1 /RYx'?± Sx i± 2na ' + rf + (n/R) 4 x 

(40) 

In the above equation The factor (1 — v) comes from \D(v, 9 — 7r)|. 

On the other hand, tension of F-string is derived from the following action, 

Uf = £™dxe*' 2 y/rl + (r/R)*. (41) 
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where r x = d r x. Then the tension of the F-string is obtained as, 

5U F 1 e^/VW 



X 



5t x 2na' ^(W + ^ i/R y 



(42) 



To compare tensions, we take a "vertical limit", namely, r[ — > oo, — > oo. Then 
the following equality holds; 

The tensions in the ^-direction vanish in the vertical limit. The above equality means 
the tension of the cusp equals to (1 — v)N times that of F-string automatically in the 
vertical limit in the case of q > 0. 

Other case except vertical limit, the following no-force condtion, 

SU „ SAT SU F 5U „ , AT 5U F 

iFr {1 ' v)N ^' terr (1 - v,N s^' (44) 

assures the relation between tensions of cusps and F-strings. 
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